function [C]=ShannonCost3s(mu,p,gamma_g,gamma_b)

C=-1*(H(mu(1))+ H(mu(2))+H(mu(3))- p*(H(gamma_g(1))+H(gamma_g(2))+H(gamma_g(3)))-(1-p)*(H(gamma_b(1))+H(gamma_b(2))+H(gamma_b(3))));
end

